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Abstract. We have discovered a new AM Her system 
as the optical counterpart of the ROSAT All-Sky-Survey 
source RXJ1016. 9-4103 (= 1RXS J101659.4-410332). The 
X-ray spectrum is very soft and the X-ray intensity is 
strongly modulated with the orbital period. Optical pho- 
tometric and spectroscopic follow-up observations reveal 
a synchronously rotating binary with an orbital period 
of 134 min, placing RX J1016. 9-4103 in the period gap. 
The strength of the TiO bands suggests a secondary spec- 
tral type later than M3 V and a distance of 615±150 pc. 
Based on two clearly visible broad humps in the optical 
spectrum (interpreted as cyclotron features) a magnetic 
field strength of 52 MG is deduced thus proving the polar 
classification. 

Key words: X-rays: stars - stars: cataclysmic variables - 
accretion, accretion discs - stars: magnetic fields - stars: 
individual: RXJ1016. 9-4103 = 1RXS J101659.4-410332 



1. Introduction 

AM Her type variables are a subgroup of cataclysmic vari- 
ables (CVs) in which the magnetic field of the white dwarf 
controls the geometry of the material flow between the 
main-sequence donor and the white dwarf primary as well 
as synchronizes the white dwarf spin period with the bi- 
nary orbital period (see e.g. Warner 1995 for a detailed re- 
view) . The inflow of matter along the magnetic field lines 
(of one or occasionally also two magnetic poles) is decel- 
erated above the white dwarf surface producing a shock 
front. This region is thought to emit hard X-rays (usu- 
ally modelled in terms of thermal bremsstrahlung of 10- 
20 keV) and polarized cyclotron radiation (hence these 
systems are named also polars) in the IR to UV range. 
In addition, a strong soft component has been frequently 
observed from polars which is thought to arise from the 
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* Partly based on observations collected during MPI time at 
the 2.2 m telescope of the European Southern Observatory, La 



Table 1. Log of observations 



Telescope" 



Date 



Energy Band 



r int 
(sec) 



X-ray 

ROSAT P 1990 Dec 09-11 0.1-2.4 keV 

ROSAT P 1993 Nov 30-Dec 4 0.1-2.4 keV 

ROSAT P 1994 Jun 01 0.1-2.4 keV 

ROSAT H 1997 Dec 08 0.1-2.4 keV 

optical spectroscopy 

ESO 2.2 m 1995 Mar 26 3000-9000 A 

optical photometry 

ESO 2.2 m 1995 Mar 25-27 B, R 



430 
6760 
2930 
1800 

1800 

600, 180 



a The abbreviations have the following meanings: ROSAT P 
or H = PSPC or HRI detectors onboard the ROSAT satellite, 
ESO = European Southern Observatory, La Silla/Chile. 



heated accretion pole (usually modelled in terms of black- 
body emission with 20-50 eV temperature). 

It is this soft X-ray component which has led to the 
discovery of a few dozen new polars by ROSAT obser- 
vations over the last 8 years, most notably the ROSAT 
all-sky survey (e.g. Beuermann & Burwitz 1995 for a first 
summary). Based on these discoveries, the strength of the 
soft component (or more precisely the ratio of the black- 
body and the bremsstrahlung component) has been found 
to increase with the magnetic field strength of the white 
dwarf. Since simple reprocession of the hard component 
cannot produce the observed soft X-ray fluxes, a scenario 
of high-m, blobby accretion has been proposed as a mech- 
anism to explain the soft excess (Kujpers & Pringle 1982). 
In a different approach, bremsstrahlung emission is grad- 
ually replaced by more efficient cyclotron cooling at high 
magnetic field strengths (Woelk & Beuermann 1996). 

The source described here has been discovered as a re- 
sult of a systematic survey for supersoft X-ray sources in 
the all-sky survey data (see Greiner 1996 for details of this 
survey) which revealed a large number of CVs and single 
white dwarfs. Other polars identified from this sample in- 
clude V844 Her = RX J1802. 1+1804 (Greiner, Remillard 
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Fig. 1. A 10 min 5-band CCD image of RXJ1016. 9-4103 
obtained with the 2.2 m telescope at ESO. The circle 
denotes the 3cr X-ray error circle derived from the HRI 
pointing (10" radius). The optical position is measured as: 
R.A. (2000.0) = 10 h 16 m 58 i ;9, Decl. (2000.0) = -41°03'45" 
(±1"). The bright star at the left denoted with "C" has 
been used as comparison star for the photometry. 

witz etal. 1996), and RX J1724. 0+4114 (Greiner, Schwarz 
and Wenzel 1998). In this paper we present photomet- 
ric, spectroscopic and X-ray observations (summarized 
in Tab. [j]) which led to the discovery of another polar, 
RXJ1016. 9-4103 (henceforth referred to as RX J1016). 

2. X-ray discovery and identification 

RX J1016 was scanned during the ROSAT all-sky-survey 
over a period of 4 days in December 1990 for a total ob- 
serving time of 430 sec. Its mean count rate in the ROSAT 
position-sensitive proportional counter (PSPC) was 0.13 
cts/s, and the hardness ratio HRI = —0.95 ± 0.12 where 
HRI is defined as (H-S)/(H+S), with H (S) being the 
counts above (below) 0.4 keV over the full PSPC range of 
0.1-2.4 keV. 

For the timing analysis the source photons were ex- 
tracted with a radius of 4'. The background was chosen at 
the same ecliptic longitude at «1° distance, correspond- 
ing to background photons collected typically 15 sec before 
or after the time of the source photons. Standard correc- 
tions were applied using the dedicated EXSAS software 
package (Zimmermann etal. 1994). The RASS light curve 
folded over the best-fit period as derived in section 3 is 
shown in the upper panel of Fig. ||. The X-ray flux shows 
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Fig. 2. Results of the period-search using the analysis- 
of-variance method. The adopted period of 134±3 min is 
marked as the shaded region. 



and a pronounced long faint-phase where the X-ray flux 
is nearly zero (formal count rate of 0.05 ± 0.05 cts/s). 

The soft X-ray spectrum and strong variability imme- 
diately suggested the cataclysmic variable nature. Indeed, 
using the best-fit X-ray coordinate derived from the all- 
sky survey data only one optical object brighter than 20 m 
was found within the error circle (see Fig. |l|) which later 
was identified spectroscopically as a polar (see section 4). 

Dedicated follow-up pointed ROSAT observations 
were performed in November 30 - December 4, 1993 and 
June 1, 1994 with the PSPC and on December 8, 1997 with 
the high-resolution imager (HRI) . Table ^| summarizes the 
relevant numbers of these measurements including the to- 
tal number of detected counts (column 3), the mean count 
rate (column 4), the hardness ratio (column 5), the cov- 
erage of the orbital phases (column 6), the off-axis angle 
(column 7) and the best-fit X-ray positions (column 8). 
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Table 2. Basic numbers of the X-ray observations of RXJ1016. 9-4103 



Date 




Obs-ID (1) 


N c t s 


mean count 


HRI 


phase 


offaxis 


X-ray position 












rate (cts/s) 




coverage 


angle 


(equinox 2000.0) 




1990 Dec 09- 


-12™ 




57 


0.13 


-0.95±0.12 


5% 


0'-52' 


10 h 16 m 58:2-41 o 03 / 38" 


± 15" 


1993 Nov 30- 


-Dec 4 


201605p 


430 


0.06 


-0.98±0.01 


55% 


0^50 


10 h 16 m 59 ! :8-41 o 03 , 52" 


± 25" 


1994 Jun 1 




201605p-l 


352 


0.12 


-0.98±0.01 


30% 


0^32 


10 h 16 m 59 ! :2-41 o 03'44" 


± 25" 


1997 Dec 08 




300580h 


3 


0.0019 (3) 




20% 


0'.21 


10 h 16 m 58 , :9-41 o 03'40" 


± 10" 



' ' The letters after the number denote the detector: h = HRI, p = PSPC. 

( 2 ' For the position determination only photons in the energy range 0.25-0.5 keV have been used to avoid position degradation 
due to ghost images. 

^ Note the lower sensitivity of the HRI at soft energies by a factor of 7.8 as compared to the PSPC. 



3. Optical and X-ray variations 

CCD photometry was obtained during 3 nights in March 
1995 with EFOSC II at the 2.2 m telescope on La Silla 
(during MPI time). Observations were performed mostly 
with the B filter, though a few exposures also were made 
with the V and R filters. Exposure times ranged between 
10-600 sec. Observational details are listed in Tab. § The 
images were processed using the profile-fitting scheme of 
the DoPhot reduction package (Mateo & Schechter 1989) 
to achieve high accuracy. The optical light curve is charac- 
terized by a broad bright phase and a shallow 0.3 mag deep 
faint phase. In contrast, the X-ray intensity variations are 
much more pronounced and show a clear on/off- like be- 
haviour. 

In order to derive an orbital period we carried out 
a period search using the analysis-of-variance method 
(Schwarzenberg-Czerny 1989) for both, the optical pho- 
tometry data as well as the combined ROSAT PSPC 
pointings and the phase-dispersion method (Stellingwerf 
1978) for the ROSAT all-sky survey data. The resulting 
periodograms (Fig. ^) show a variety of maxima which are 
predominantly caused by the poor sampling rather than 
the intrinsic variability. A more closer look at the bottom 
two panels (optical and ROSAT pointing data) suggests 
that common periods of 123, 134 or 149 min are possi- 
ble (134 min is the most reliable period from the ROSAT 
pointed data alone) . Inclusion of the ROSAT all-sky sur- 
vey data clearly favours the 134 min period though we cau- 
tion that due to the orbital motion of the ROSAT satellite 
(96 min) the signal may be affected. We note also that the 
actual best-fit period of the -ROSAT all-sky survey data 
is 132.5 min, thus being slightly smaller than the period 
derived from the optical and pointed ROSAT data. Note 
that for each data set (corresponding to the three panels 
in Fig. ||) a slightly different best-fit period is obtained 
though the errors (which we consider to correspond to the 
width of the peak) overlap. Overall, based on the appear- 
ance of the folded light curves (Fig. ||) according to a cri- 
terion of "greatest simplicity" we feel that P or b = 134±3 
min is our best estimate (shaded area in Fig. ||). 



4. A low resolution spectrum 

4.1. Optical identification 

The identification spectrum of the likely counterpart was 
obtained on March 26, 1995 with the ESO 2.2 m telescope 
at La Silla/Chile and was exposed for 30 min. We used the 
EFOSC II spectrograph equipped with a lkxlk Thomson 
CCD detector with a 450 A/mm grism (corresponding to 
about 8.5 A/pixel) covering the optical wavelength range 
from 3000-9000 A. The observation was performed under 
stable photometric conditions and accompanied by mea- 
surements of the standard star LT4816. We additionally 
applied a correction using the B magnitude derived from 
a direct image taken just prior to the spectroscopic obser- 
vation to calibrate the flux with an accuracy of ~ 10% 
(using standard MIDAS procedures). By convolving the 
original spectrum with functions representing the BVR 
bandpasses we arrive at B = 18™3, V = 18™5 mag and 
R = 18'. n with a mean error of ±0.2 mag for RX J1016. 

The resulting spectrum plotted in the upper panel of 
Fig. H reveals strong emission lines with of the Balmer 
series and He 11 typical for a CV. The high-excitation line 
Hen A4686 A is quite prominent. Its strength (2/3 of the 
H/3-flux) and the inverted decrement of the Balmer lines 
clearly indicate a magnetic nature of the object. 

4.2. The secondary star 

We further deconvolved the optical spectrum by fitting a 
power-law in the blue range between 3800-5000 A. After 
subtracting this blue continuum contribution the late-type 
companion is revealed by the slope of the continuum in the 
near-IR and obvious flux depression due to the TiO ab- 
sorption throughs which are best seen in the AA 7053 and 
7667 A bandheads. The faintness of the features prevents 
an accurate quantitative estimate of the spectral type, but 
the low ratio Si 0.5 of the flux deficits Fxio(6165A) to 
F T io(7667A) suggest that it is later than dM 3. The ex- 
pected parameters for a Roche-lobe filling secondary in a 
134 minute system are M 2 = 0.18M Q , R 2 = 0.22R Q and 
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Fig. 3. Optical and X-ray light curves of RX J1016 plot- 
ted over orbital phase for a period of P = 134.4 min. 
Note that the phase designation of the middle two panels 
is arbitrary with respect to the all-sky survey (top panel) 
and optical (bottom panel) data. Data are plotted twice 
for clarity and units are PSPC cts/s for the upper three 
panels and B magnitudes for the bottom panel. 

mass-radius and mass-spectral type relations of Patterson 
(1984) and Kirckpatrick & McCarthy (1994) were used. 

The distance of RX J1016 can be estimated from the 
observed strength of the M-dwarf using the method of 
Bailey (1981). In the following we take the dM 4.5 star 
Gl 83.1 (K = 6™67) as a template and assume that both 
M-stars have the same colour V — K = 5™62 (Leggett 
1992). By comparing the measured flux deficits in the 
A7667-band in Gl 83.1 and that of RXJ1016 (1800 10 16 
and 0.2 10 16 erg cm~ 2 s _1 ), we derive a K magnitude of 
16™5 for the secondary in RX J1016 . If we take this value 
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Fig. 4. upper panel: Low resolution optical spectrum 
of RXJ1016.9-4103 obtained on March 26, 1995. The 
dashed line represents a power-law with f\ oc A~ 2 fitted in 
the range between 3800-5000 A. middle panel: Residual 
spectrum after subtraction of the blue continuum shown 
together with the spectrum of the late M-dwarf Gliese 
83.1 (M 4.5) scaled to match the strength of the observed 
TiO features, lower panel: The spectrum after the re- 
moval of the contribution of the secondary star. It exhibits 
broad emission features at 5500 and 7400 A, tentatively 
identified as 3rd and 4th harmonic of the cyclotron fun- 
damental. The solid line represents our favoured modcll 
with B = 52 MG, 9 = 80° and for kT = 10 kcV. The 
alternative interpretation for an B = 41 MG field is also 
given (dashed line). 

Beuermann & Weichhold (1998) a distance of 615 ± 150 
pc is deduced. The error includes an uncertainty of 25% 
for both, the flux calibration and the scaling of the TiO 
band deficits. 

Detailed phase-resolved spectroscopic studies of se- 
lected systems, e.g. AM Her (Davey & Smith 1996), QQ 
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Fig. 5. Left: Phase-averaged X-ray spectrum of the two merged PSPC observations of RX J1016 unfolded with the 
sum of a blackbody and a thermal bremsstrahlung spectrum with a temperature fixed at 20 keV (see text for details). 
The lower left panel shows the residua of the fit in units of a. Right: The 68%, 90% and 95% significance contours of 
the blackbody model fit to the merged PSPC spectrum of RX J1016 in the Nn — kT plane. Overplotted are the lines 
of constant luminosity in units of (D /100 pc) 2 (dashed), as well as the value of the total galactic absorbing column 
(vertical dotted line) according to Dickey & Lockman (1990). 



etal. 1998) have demonstrated that absorption features 
from the illuminated side of the secondaries in magnetic 
CVs can be heavily suppressed due to the strong X/UV 
irradiation from the accretion region. We therefore cau- 
tion that the distance derived from our single spectrum 
with unknown viewing geometry of the binary system can 
serve as an upper limit only. 

4.3. A magnetic Geld estimate 

After subtraction of a suitably scaled spectrum of Gl 83.1 
two broad emission humps at AA 5500 and 7400 A clearly 
stand out in the red part of the spectrum (Fig. [|, lower 
panel), which we tentatively identify as cyclotron lines 
from the hot optically thin post-shock region. 

An unambiguous estimate of the magnetic field 
strength is difficult to achieve with only two distinct 
humps observed. Their separation favours an interpreta- 
tion as the 3rd and 4th harmonic of the cyclotron fun- 
damental, in which case the field strength is 52±2 MG 
slightly depending on the plasma temperature and the yet 
unknown field orientation. An alternative interpretation 
as the 4th and 5th harmonics, in which case B is ~ 41 MG, 
cannot not completely excluded but results in a distinct 
mismatch of the expected humps with the shape of the 
continuum between 6000-7000 A. Normalized cyclotron 
absorption coefficients for both cases assuming 9 = 80° 
and kT = 10 keV shown in Fig. [| illustrate this ambigu- 
ity. 



5. The X-ray spectrum 

The two pointed PSPC observations have resulted in the 
detection of nearly 800 photons thus enabling spectral in- 
vestigation. For the spectral analysis the source photons 
of the two PSPC pointings were merged and then ex- 
tracted with a radius of 1.'5. The background was chosen 
from concentric circles around the source region with a 
radius of 3'. Other nearby sources were cut out, and the 
background area normalized to the source extraction area 
before the background subtraction. Standard corrections 
were applied using the dedicated EXSAS software package 
(Zimmermann etal. 1994). As a first step, we have con- 
sidered the total spectrum as collected during each of the 
two observations (after merging). Adopting a blackbody 
plus a thermal bremsstrahlung model (the temperature of 
the latter fixed to 20 keV) we find that there is no need 
at all for the inclusion of the bremsstrahlung component, 
that is the spectral fit does not require any spectral com- 
ponent for hard emission (Fig. [s]). The resulting best-fit 
temperature of the blackbody component is 10±7 eV and 
the absorbing column ]Vh=6.4x 10 20 cm -2 , e.g. about 50% 
of the total galactic absorbing column in this direction 
(Dickey & Lockman 1990). However, this gives an unreal- 
istic high luminosity for a cataclysmic variable due to the 
very low temperature. If we fix the temperature at the 
canonical value of 20 eV, the best-fit absorbing column is 
only ]Vh=2.4x10 20 cm" 2 . The difference in reduced \ 2 1S 
only marginal for all temperatures in the 10-25 eV range. 
As a second step we have investigated the emission dur- 
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for a total of 7340 sec. Selecting time intervals where the 
60 sec averaged X-ray intensity is less than 0.05 cts/s, we 
find that the residual emission is consistent with being 
background radiation. The mean count rate during the 
X-ray bright phase (at >0.05 cts/s) is 0.38 cts/s in the 
PSPC. 

Considering only the X-ray bright phase, the unab- 
sorbed fluxes of the two model components in the ROSAT 
band (0.1-2.4 keV) arc F hhdy = 5.1 xlO" 11 erg cm" 2 s" 1 
(using kTbbdy=20 eV for comparison purposes) and -Fthbr 
<9.8xl0 -13 erg cm~ 2 s _1 , giving a small flux ratio of 
-Fthbr/-Pbbdy < 0.02. The mean bolometric luminosity dur- 
ing the X-ray bright phase is (again with kTbbdy=20 eV) 
L x = 2.4xl0 32 (D / 100 pc) 2 erg/s. 

6. Discussion and Conclusion 

The measured period of P = 134 min places RXJ1016 
in the 2-3 hr CV period gap which is thought to be due 
to the transition from orbital angular momentum loss by 
magnetic braking (P> 3 hr) to gravitational radiation 
(P< 2 hr) (King 1988). Over the last years, more polars 
have been found with periods in the 2-3 hr range thus 
re-initiating the debate on the existence and significance 
of the period gap for magnetic systems (Wickramasinghe 
& Wu 1994, Wheatley 1995). 

The shape of the X-ray light curve strongly argues for a 
self-eclipsing polar. The length of the X-ray bright phase of 
^0.5 phase units suggests (within our limited accuracy) 
an one pole accreting geometry such that the accretion 
region passes behind the limb of the white dwarf for half 
the orbital period. The lack of eclipses implies i < 78°. 

Our tentative estimate of the magnetic field strength 
is 52±2 MG and the strong soft X-ray excess are in good 
agreement with the correlation found between these quan- 
tities (Beuermann 1998). 

The strongly modulated soft X-ray emission, the seem- 
ingly synchronous rotation, the strength and relative in- 
tensities of the Balmer and He n emission lines as well as 
the cyclotron humps in the optical spectrum clearly sug- 
gest the polar nature of RX J1016 though no polarimet- 
ric measurements have been made. More detailed optical 
photometry and phase-resolved spectroscopy as well as po- 
larimetry are needed to determine the system parameters 
of this new polar. 
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